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© An MRI apparatus of the invention comprises a magnet system for producing a homogeneous magnetic field 
along the z axis. A plurality of wire conductors are arranged in a "bird cage" or "squirrel cage" configuration and 
are supplied with a circularly rotating RF field centered near the Larmor frequency o> 0 . The wire conductors also 
receive a voltage source V p having a frequency o p , which is less than the Larmor frequency. The voltage source 
V p is rotated about the wire conductors at a frequency of o s , which is much less than the frequency *> p , to 
successively apply the voltage V p to each pair of diametrically opposed wire conductors in a manner which is 
ordered as a spatial cylindrical harmonic rotating at a frequency « s . The voltage source V p produces electric 
fields in the object which, in turn, generate linear magnetic field gradients in the object. The voltage source V p 
phase encodes the signals detected from the object which thereby enables a transverse magnetization at each 
point within the object to be rapidly mapped back by projection techniques. Alternately, the transverse 
magnetization can be mapped back by an inverse transformation of a complex plane extending from the 
measured region by analytic continuation or by a linear transformation utilizing the solutions to So mm erf eld's 
Integral. Also, the complex permittivity or ionic mobility may be continuously measured through its effect on the 
transverse magnetization for analysis or process control. With the invention, an image of the object is rapidly and 
continuously produced without requiring the use of the conventional gradient field coils and without any 
mechanical moving parts. 
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This invention is generally concerned with magnetic resonance imaging (MRI) and, more particularly, to 
the formation of gradient fields in an imaging volume. 

A prior art MRI apparatus, such as that depicted in Fig. 1 and disclosed in U.S. patent No. 4,952,878, 
which is hereby incorporated by reference, is comprised of a magnet system 2 and a magnet supply 
5 source 6 for generating a homogeneous stationary magnetic field H D . A magnet system 4 and a supply 
source 8 generate magnetic gradient fields and an RF coil 10, connected to an RF amplifier 12, generates a 
radio frequency alternating magnetic field. 

The uniform magnetic field H G produced by the magnet system 2 and the magnet supply source 6 is 
formed within an imaging volume along a z axis and causes the nuclei within the object to tend to be 
w aligned along the z axis. The nuclei precess about the z axis at their Larmor frequencies according to the 
following equation: 

w 0 = y H 0 EQ. (1), 

75 where o> 0 is the Larmor frequency at a field strength of H 0 and y is the gyromagnetic ratio. The 
gyromagnetic ratio y is a constant and is a property of the particular subatomic particle. For instance, the 
gyromagnetic ratio y for water protons is 4.26 kHz/gauss. Thus, in a 1.5 Telsa polarizing magnetic field H Q , 
the resonant or Larmor frequency of water protons is approximately 63.9 HHz. 

In a typical imaging sequence, the RF coils 10 are used to generate an RF signal centered at the 

20 Larmor frequency o> 0 . In conjunction with this RF signal, a magnetic field gradient G 2 is generated with the 
gradient coils 4 so that only the nuclei in a selected slice through the object along the x-y plane, which have 
a resonant frequency centered at the Larmor frequency w Q , are excited and absorb energy. Next, the 
gradient coils 4 are used to produce magnetic field gradients along the x and y axes. Generally speaking, 
the x axis gradient G x spatially encodes the precessing nuclei continuously by frequency and the y axis 

25 gradient G y is used to pulse phase encode the nuclei. 

The RF coils 10 may be formed in the design of a "bird-cage" or a "squirrel cage" as shown in Fig. 2. 
Other RF coil designs are also known, such as those disclosed in U.S. Patent No. 4,717,881, U.S. Patent 
No. 4,757,290, U.S. Patent No. 4,833,409, U.S. Patent No. 4,916,418, U.S. Patent No. 5,017,872, U.S. Patent 
No. 5,041,790, U.S. Patent No. 5,057,778, and U.S. Patent No. 5,081,418. The RF coils 10 generate a 

30 circularly polarized RF field which rotates near the Larmor frequency over the antenna wires 30 of the coil. 
The RF coils 10 are also used to detect signals from the object and to pass the detected signals onto a 
signal amplifier 14, a phase-sensitive rectifier 16, and to a control device 18. 

A description of_the above process is given in more detail_below. The spins of the nuclei have an 
angular momentum S which is related to their magnetic moment u through the gyromagnetic ratio 7, by the 

35 equation: 

U = y S EQ. (2). 

With reference to Fig. 3, a single nucleus 1 has a magnetic moment 3. All of the spins are subjected to a 
40 total ambient magnetic field H defined by the equation: 

H = lh x + ~jh y + k(H 0 + h 2 ) EQ. (3). 

As a result of the total magnetic field H, a torque T is produced on each spin. With reference to Fig. 3, the 
45 total magnetic field H produces a torque T on the nucleus 1 to cause the nuclei 1 to precess. The torque T 
produced is defined by the following equation: 



T = u> xH = 

50 




EQ. (4). 

55 Substituting and rearranging from the above equations yields: 



2 
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EQ. (5), 



or, 

10 di = (* x y H)dt EQ. (6). 

From equation (6), it is evident that the differential spatial phase increment di> of each spin is related to a 
cross product of the direction of each spin <l> with the ambient magnetic field H multiplied by the 
gyromagnetic constant -y. 
75 Substituting from the above equations yields: 

di = (¥ x y (ih x + )h y ))dt + (¥ x y k(H 0 + h z ))dt EQ. (7), 

or 

20 



cflF & dB + cfij> 

25 EQ. (8), 

where de is the differential increment of spatial phase of the precession of u about an axis perpendicular to 
the z axis and d<f> is the differential increment of the spatial phase of the precession of u about the z axis 7. 
30 Therefore, the time-varying angular relationship of each magnetic moment u is controlled by the applied 
ambient magnetic field H. 

A differential increment 60 of u is produced when ih x + jh y rotates synchronously with $ at the Larmor 
frequency. Often, the RF coils 10 are used to create such a circularly polarized RF magnetic field which 
rotates at the Larmor frequency « 0 by selectively choosing the values of ih x + jh y . As a result, an 
35 interaction energy W m is transferred to the spins and is given by: 

dW m = -T*d$ = -(Uxk(H 0 + h z )).d$ EQ. (9), 

W m = u(H 0 + h 2 )(cose - 1) EQ. (10). 

40 

If the spins are coupled to the surrounding medium to permit relaxation to the minimum energy W m 
state before the application of ih x + jh y , then the average value erf 0 for a discrete volume of spins will be at 
a minimum. The vector sum of each of the magnetic moments u in a discrete volume will then produce a 
net magnetization M aligned with the z axis. The precessions of the spins d produced by ih x + jh y about an 
45 axis perpendicular to the z axis will give the net magnetization M a transverse component M T with a spatial 
phase <f>. The strength of the transverse magnetization M T is a function of the number of spins in each 
discrete sample^ volume and of the previous excitation and relaxation history of these spins. The transverse 
magnetization M T is a component of the intensity of a total rotating flux density B T through the permeability 
of free space u D where: 

50 

B T = u 0 0h x + 'jhy + Mr) EQ. (11). 

If an antenna loop is placed with its area vector perpendicular to the z axis, it will have a voltage induced in 
it which is proportional to the time derivative of B T and also a function of the geometry of the conducting 
55 antenna loop and of the spatial distribution of the discrete sample volumes, each containing a singular value 
of B T . This functional dependance is given by the Maxwell-Faraday equation and the vector reception field 
of the antenna geometry, which may be determined by experimentation. 

The transverse magnetization M T has an instantaneous rotational frequency: 
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75 



40 



45 



55 



H = y(H 0 * h 2 ) a o> 0 + oi z 



EQ. (12) 



10 



EQ. (13). 

equ^tio P n° lar ***** fie ' d @ " defin8d 38 ^ 9r3dient ° f hz 3 " d ' S determined according to the following 



20 G = V£ = l — i + j _£ + Tc -=-£ 

dx dy 3z 



EQ. (14). 



25 As discussed above with reference to Fig. 1, the gradient magnets 4, or other type of "gradient coils - 
prov.de a spatial address for each discrete component of the imaging volume. The polar vector field G" 
having components G x , G y , and G 2 defined by Equation 14, is successively incremented by driving the 
gradient coils over a longer pulse interval or by incrementing the amplitude of the driving current The 
signal from the imaging volume is analyzed by Fourier transformation or other integral transformation to 

30 produce a mathematical representation of the spatial address of the transverse magnetization M T in the 
complex plane. 

A problem with many prior art MRI systems, however, is that the time needed to produce the image is 
relat.vely long, for instance about 6 minutes. During this time the object to be imaged, for instance the 
patient, must remain motionless. " ' 

as The prior art MRI systems are also at a disadvantage due to the manner in which gradient fields are 
produced Gradient field generating systems, such as the supply system 8 and the gradient magnet system 
4 of Fig. 1. must be able to quickly pulse large currents while precisely controlling the magnitudes of these 
currents. The pulsing of the large currents produces a pounding noise which often times must be blocked 
out with earplugs. Although many types of MRI pulse sequences can help reduce the time needed to scan 
an .mage, the circuitry in these supply systems becomes extremely complex as well as very costly in order 
to provide precise control of large currents. 

r J*L additi0na L pr ° blem With many MRI sysfems is that the V are ""able to generate rotating gradient 
fields that are sufficiently linear for projection reconstruction image acquisition techniques Thus these MRI 
systems are incapable of producing images by techniques such as projection reconstruction, but must 
instead rely upon a two or three dimensional truncated Fourier series and the Fourier transformation to 
produce an Image. 

The prior art MRI systems further suffer from a disadvantage in that the gradient fields are only 
produced to a first order approximation of the desired linear gradient fields. Because the gradient fields 
produced only approximate the desired linear fields, some of the image reconstruction systems used in the 
so prior art systems suffered as a result. 

Accordingly, it is an object of the invention to provide an MRI apparatus which overcomes these as well 
as other problems in the prior art. 

It is an object of the present invention to provide an apparatus which can produce linear magnetic field 
gradients. 27 

It is also an object of the present invention to provide an apparatus that produces magnetic field 
gradients by generating time varying electric fields. 

It is a further object of the present invention to provide an apparatus which generates a "real-time- 
image of an object. 
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It is another object of the present invention to provide an apparatus which generates an image of an 
object by projection reconstruction techniques. 

It is yet another object of the present invention to provide an MRI apparatus that can generate an image 
without having to rely upon the two or three dimensional Fourier transformation but which uses alternate 
5 methods of processing signals detected from an object. 

It is yet a further object of the present invention to provide an apparatus for generating magnetic field 
gradients which is useful in spectroscopy as well as in other applications involving magnetic resonance. 

To achieve the foregoing and other objects, in accordance with the present invention, as embodied and 
broadly described herein, the apparatus of the present invention comprises a magnet system or other 
w conventional system for generating a homogeneous magnetic field directed along the z axis. A plurality of 
wires are arranged parallel to each other along the z axis in a typical "bird cage" or "squirrel cage" 
configuration which is tuned near the Larmor frequency c> 0 . 

A circularly polarized RF field centered near the Larmor frequency <u 0 is generated so that it rotates 
about the wire conductors. A sinusoidal voltage source having a frequency <d p , which is less than the Larmor 
75 frequency o> 0 . is sequentially applied as a spatial cylindrical harmonic across diametrically opposed wire 
conductors at an angular spatial frequency of <* s , where <d s is very much less than the frequencies a> 0 and 

The electric current density that is generated by the voltage source of frequency o> p produces a 
magnetic field in the object^ that alters the phase of each spin in a periodic and spatially ordered field. The 

20 transverse magnetization M T of the nuclei produces an RF signal that is detected by the wire conductors 
and is then supplied to a control device for signal processing. 

The control device can analyze the signal as a high signal to noise rapidly repeated continuous analytic 
function at the origin of the complex "K" plane by analytic continuation. Also, the control device may 
analyze the received signals as a continuous set of rapidly rotating views of the sum of the transverse 

25 magnetization M T by projection reconstruction. The determination of the transverse magnetization M T may 
be performed with Bessel functions or with the solution to Sommerfeld's integral, both of which allow the 
use of look up tables. 

The invention can therefore generate the necessary magnetic field gradients by applying time varying 
electric fields to an RF coil assembly. Thus, the conventional "gradient coils" are no longer a necessity in 

30 the formation of images with an MRI apparatus. However, the apparatus of the present invention does not 
necessarily eliminate the conventional gradient coils since the application of the time varying electric fields 
may be used in conjunction with the conventional gradient coils. For instance, the gradient fields produced 
by the electric fields may be used to dephase or enhance some of the phase encoding generated with the 
conventional gradient coils. As a result of this, the image generated by the MRI apparatus can be made 

35 more narrow and thus the apparatus can focus in on a particular region of interest and block out all other 
regions. 

Additional objects, advantages and novel features of the invention will be set forth in the description 
which follows, and will become apparent to those skilled in the art upon reading this description or 
practicing the invention. The objects and advantages of the invention may be realized and attained by the 
40 appended claims. 

The accompanying drawings, which are incorporated in, and form a part of, the specification, illustrate 
some embodiments of the present invention and, together with the description, serve to explain the 
principles of the invention. In the drawings: 

Fig. 1 is a block diagram of an embodiment of a prior art MRI apparatus; 
45 Fig. 2 is a prior art RF coil assembly; 

Fig. 3 is a perspective representation of a precessing nuclear spin; 
Fig. 4 is a block diagram of a first embodiment of the invention; 

Fig. 5 is an exemplary diagram of the circuitry used for generating the magnetic field gradients by way 
of electric fields; and 

50 Fig. 6 is a second embodiment of the MRI apparatus of the present invention. 

Reference will now be made in detail to the preferred embodiments of the invention, examples of which 
are illustrated in the accompanying drawings. 

With reference to Fig. 4, an MRI apparatus of the invention comprises a magnet system 2 for 
generating a homogeneous stationary magnetic field H D and a supply source 6 for supplying the magnet 
55 system 2 with operating power. The magnet system 2 may be of the superconducting type which requires a 
connection to a cooling device 26 through a cooling duct 27. The invention, however, is not limited to this 
specific type of magnet system 2 and other types of magnet systems 2, such as permanent magnets or 
resistive nonpermanent magnets, may alternatively be used. 

5 
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A phase-sensitive rectifier 16 receives a signal from an oscillator 24 and supplies an output to a control 
device 18. The control device 18, together with a modulator 20 which is under the control of the control 
device 18, produces a signal modulated near the Larmor frequency which is amplified with an RF amplifier 
12 and then supplied to an interface circuit 15. 

The interface circuit 15 may comprise an amplitude and phase control network 42 such as that 
illustrated in Fig. 2. The interface circuit 15 individually controls the voltages applied to each of a plurality of 
wires 30 forming part of an RF coil assembly 13. Preferably, the interface circuit 15 generates a circularly 
polarized RF f.eld having a frequency centered near the Larmor frequency «„ and which rotates over the 
wires 30 of the RF coil assembly 13. 

The control device 18, together with a second modulator 17 controlled by the control device 18 also 
produces a sinusoidal voltage at a frequency of «, p . This signal is amplified by a low frequency amplifier 19 
before being supplied to the interface circuit 15. The interface circuit 15 comprises circuitry 50 shown in 
diagrammatical form in Fig. 5, for generating electric fields from the sinusoidal voltage. The circuitry 50 
comprises a polyphase synchro generator 52 for receiving across its terminals the voltage source V D which 
is preferably a sinusoidal voltage signal equal to V p sin<«„t). The polyphase synchro generator 52 rotates at 
a frequency of «. and successively applies the voltage signal V p to each pair of opposed conductor wires 
30 as a spatial cylindrical harmonic function of the rotating frequency «, of the primary winding illustrated in 
rig. 5. 

The coil assembly 13 having the opposed wires 30 may be formed in the conventional "bird cage" or 
squirrel cage" configuration, such as that illustrated in Fig. 2. The coil assembly is tuned near the Larmor 
frequency « 0 so that the circularly polarized RF field generated near the Larmor frequency « c can rotate 
about the wires 30. A first set of isolation filters Z, are for filtering out the circularly polarized voltages 
generated near the Larmor frequency <o 0 and are connected between the wires 30 and contact points for the 
-Poly£!32Se_synch I og^rator 52. Another set of filters Z 3 are connected between the wires 30 on one rinq 
conductor 34 of thTcBTas^embly 13. These filters Z 3 tune the coil assembly 13 near the Larmor frequency 
*c , so that voltages near the Larmor frequency may rotate about the coil assembly 13 and also prevent 
other signals, such as the voltage signal V p , from rotating about the coil assembly 13. 

The Larmor frequency «„, the frequency <o p of the voltage signal V p , and the frequency « s at which the 
polyphase synchro generator rotates, are selected such that: 

u>„>o> p »« s EQ. (15). 

The electromechanical polyphase synchro generator configuration 52 is diagrammatic only and is 
preferably comprised of solid state circuitry. For instance, the polyphase synchro generator 52 may be 
replaced with an amplitude and phase control network and a voltage divider network that produces a spatial 
cylindrical harmonic array of sinusoidal voltages having a temporal frequency of « p .The voltage divider 
network may then periodically apply the array of sinusoidal voltages to opposed wire conductors 30 at the 
angular spatial frequency of a*. These voltages are periodically applied over a spatial angle «, where: 

40 a = o s t EQ. (16). 

Besides using the coil assembly 13 to transmit electromagnetic fields that have frequencies of « and 
-p. the coil assembly 13 is also preferably used as a receiver to detect RF signals emitted by the nuclei in 
an object to be imaged. The signals detected with the RF coil assembly 13 are passed through a set of 

45 filters Z 2 for filtering out unwanted frequencies and unwanted signals, such as those intended for 
transmission or phase modulation. The detected signals are amplified with an RF amplifier 14 rectified with 
rectifier 16, and then processed in the control device 18 in order to construct an image of the object which 
is displayed on a monitor 22. 

The embodiment illustrated in Figs. 4 and 5 does not contain the conventional "gradient coils " such as 

so the gradient coils 4 in Fig. 1. Instead, the embodiment employs an RF coil assembly 13 for transmitting an 
RF pulse and for also creating the magnetic field gradients. The manner in which these gradients are 
produced and the manner in which the detected signals are processed will now be presented in more detail 
below. 

As shown by Equation 17, the gradient field G is defined to be the gradient of the z component h 2 of 
55 the magnetic field H. 



30 



35 
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EQ. (17) 

5 

As a result, with reference below to Equations 18 and 19, the curl of the gradient field G must be equal 
to zero and, using Stokes theorem, the closed line integral of the gradient field G must also be equal to 
zero. Although the curl of the gradient field G is equal to zero, the curl of the magnetic field H is not equal 
w to zero. 

VxGeO, (butVxS^O) EQ. (18); 

and 

75 

0G • dL m 0 EQ. (19); 

where, in a cartesian system of coordinates, the differential path element dL is defined as: 

20 dL = Idx + ~jdy + kdz EQ. (20). 

Therefore, the relationship between the z component h 2 of the applied magnetic field H and the gradient 
field G is as follows: 



JO JQ JO Z Ur.«) 



EQ. (21), 



and, assuming cylindrical symmetry, 

= ¥± = HHi = o 

dz dz dz 

EQ. (22). 

40 

one may therefore conclude that the relationship between the gradient field G and the magnetic field H 
is as follows: 

G = k x (V x H) EQ. (23). 

45 

The relationship between the z component h 2 of the magnetic field H and the applied electric field E 
may also be determined. From one of Maxwell's equations, it is known that the current density J is equal to 
the curl of the magnetic field H. 

50 V x H = J EQ. (24). 

As shown in Equation 25, the total current density J T is comprised of a number of components. 

7 T = ~iJ* + ]J y + kJ 2 = Ji + Jie + J7 0 EQ. (25); 

55 

where: 

Ji = impressed (source) electric current density; 
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J je - effective conductive electric current density = (o s + a> p f")E; EQ. (26) 

Jde - effective displacement electric current 
5 density; 

= j "p € ' 1; EQ. (27) 

€ = the complex permittivity e of the medium; 

w J € ' EQ. (28); 

and 

o s = the static conductivity. 

An electric field E is defined in exponential terms as follows: 

75 

E^Ee^ 

20 EQ. (29); 

where 

E = iE x +jE y EQ.(30). 

25 

Assuming a source free region, the impressed electric current density J t is zero and the total current 
density J T is reduced to the form: 

Jt = 0 + (o s + + j^)E EQ. (31). 

30 

The total current density J T expressed in more general terms is as follows: 



35 



J T a c {j „ p) e 4 \c\e*+E 

EQ. (32). 



40 By using the relationship between the gradient G and the z component h 2 of the magnetic field H (Equation 
21) and the relationship between the gradient field G and the magnetic field (Equation 23), the z component 
h 2 of the magnetic field H may be expressed in terms of the electric field E by the following equation- 



45 



so EQ. (33). 

Thus, the z component h z of the magnetic field H, which defines the gradient field G. may be generated 
by applying electric fields to the object. The relationship between^ the desired z component h z and the 
electric field E, and thus the relationship between the gradient field G and the electric field E, is determined 
55 by Equation 33. 

To aid in the understanding of this relationship, consider a uniform electric field E produced by a 
voltage V applied across two large conducting plates. The two conducting plates are perpendicular to the x 
axis, parallel to the z axis, and are separated by a distance D. The components of the electric fields 
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produced by this voltage V are then determined according to the following equations: 
E x - 3 ; E y = 0 

EQ. (34). 

io Integrating Equation 12 to solve for the phase </> as a function of time, yields the following equation: 



<t>(t> - / 0 C Y<tf o + h a )dt 



75 

EQ. (35). 

Then, substituting the value of the z component h 2 of the magnetic field H from Equation 33 and solving 
20 yields the following in time-harmonic form: 

25 

EQ. (36); 

or, in real-time form: 

30 

<t> (t) = yH Q t + y sin(o> p t * *) 

35 

EQ. (37). 

Since the transverse flux density B T produced by each element in the object having an area of dxdy in the 
complex xjy plane is known to be: 

40 

r (x,y, o r ^ r <jr.y> 

EQ. (38), 

the voltage induced in the coil assembly 13 by each element of dxdy is: 

50 

EQ. (39). 

Integrating this equation and solving for the induced voltage yields: 
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EQ. (40). 

In Equation 40, v D (t) is the output voltage detected by the coil assembly 13 which is analyzed by the 
control device 18. The element B (yit) in Equation 40 is a phase modulating function of space and time which 

w is also determined by the control device 18. Since v Q (t) and B (y t) are known variables, Equation 40 may be 
solved to determine element A^, and thus the projection of the transverse magnetization M T . 

This yields one view of a projection-reconstruction of the image. If Equation 40 is subject to a complex 
spatial transformation, such as the Fourier transformation, an analytic continuous segment in the complex 
W K" plane may be determined. 

75 Up to this point, the transverse magnetization M T has been determined as if the electric fields have 
been generated from two parallel conducting plates. Since the embodiment of Fig. 4 has a plurality of wires 
30, each pair of opposed wires 30 will approximate a pair of conducting plates for an infinitesimal element 
in the object. By performing a conformal transformation, which is within the skill of one of ordinary skill in 
the art, the generalized orthogonal rectilinear space is transformed through the Cauchy-Reimann condition 

20 equations to an orthogonal curvi-linear space so that the phase encoding electric field E may be produced 
by an array of conductors arranged parallel to the z axis around the periphery of the object to be imaged. 
The spatial address of each element M T is also subjected to the same conformal transformation. The 
purpose of the conformal transformation is to translate the relationships defined by two parallel plate 
conductors to the relationships which exist in the embodiment of Fig. 4. 

25 Next, by slowly rotating the voltage source V p about the wires 30 at the spatial angle of a, a continuous 

set of views of projection reconstruction are generated. Each view is the average of a number of images in 
time so that the signal to noise ratio is thereby improved. The rotation of the voltage source V p may also be 
visualized as generating an "n-leaved rose" pattern of closely spaced lines centered at the origin of the 
complex "K" plane. 

30 The spatial distribution of M T may be continuously computed to thereby form a virtual "real-time" image 

while simultaneously applying RF energy to or receiving RF energy from the object to be imaged. 

An alternate manner of determining the transverse magnetization M T involves the use of Bessel 
functions. Turning back to the parallel conducting plates analysis, the phase <f> as a function of time may be 
expressed by the following equation: 

35 

♦fy.t) " / 0 C Y(*o + h x )dt 

40 EQ. (41). 

After substituting in the expression for the z component h 2 of the magnetic field H, the phase function </> 
may be expressed in real time as: 

4> = yH 0 t ♦ y(y) (C) (-£) (-i) sin (w p t) 

EQ. (42). 

50 

Also, the transverse magnetic flux density B T is known to be equal to: 



EQ. (43). 



10 
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If a coil is placed in the x-z plane, a voltage v (t) will be induced in the coil which is equal to: 

EQ. (44) . 

Since Equation 42 provides an expression for <f>, sin<j> in Equation 44 may be expanded to equal: 
sin<£ = sin(o> 0 f + Ay sino p f) EQ. (45). 
In Equation 45, the variable A and o> 0 are defined to be: 



75 



A * 5^ ; * Y *° 

20 EQ. (46). 

Using Bessel functions, the expression s\n<t> may be expressed as: 



25 



sin* = J Q (Ay) sin (u> 0 t) + J^^Uy) ( -1 ) n * 1 sin ( co c * nu p ) t 



30 



(-1) "sin (g) 0 -/3(i> ) t 



35 



EQ. (47); 



40 



and 



45 



50 



•cos (o> 0 +nci) p ) t + ( -1) n (o) o -/30> p ) cos (w 0 -.no) p ) tj 



EQ. (48) . 

55 In conventional systems, a free induction decay signal is produced by broadening the absorption line 
by a conventional gradient coil. The free induction decay signal is the Fourier coefficient of the transverse 
magnetization _M T projected perpendicular to the gradient. In contrast, with the embodiment of Fig. 4, an 
electric field E modulates the phase increment to produce a line spectrum centered at the Larmor 

11 
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frequency. The amplitudes of the line spectrum are equal to the integral of the projection of M T taken along 
the electric field E weighted by a Bessel function whose order is the order of each harmonic and whose 
magnitude is proportional to the position of each element of the projection along the gradient perpendicular 
to the electric field E. 

For each harmonic n in the line spectrum, the sum of the upper and lower sideband voltage output 
detected by the coil assembly is equal to: 



*D/2 

70 ' "' J -D/2 



EQ. (49), 



15 where: 



20 



/+D/2 



EQ. (49a) 



From Equations 49 and 49a. The projection ofjhe transverse magnetization M T may be determined Once 
the projection of the transverse magnetization M T has been determined, it is within the capability of one of 
25 ordinary skill in the art to generate an image of the object. 

A third manner in which the control device 18 may process the detected signals and generate an imaqe 
of an object is through the use of Sommerfeld's integral. The transverse magnetic field density B T which is 
the output detected by the coil assembly 13, is equal to: 



30 



MO ^° J -D/2 {y) 

35 EQ. (53) . 

Taking the Fourier transform of the transverse magnetic field density B T , yields the following: 



40 



45 



50 



•e- jAyBin ^ c d(jo> p t)dy 

EQ. (54) 

If the following variables are defined accordingly: 



55 
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( O) - (a) ) 

n k - ; p a Ay; x a u> D t 

EQ. (55); 



then part of Equation 54 may be abbreviated to SommerfelcTs integral, as shown by the following: 

10 

J-d/2 °;w p J-> 

EQ. (56). 

Therefore, Equation 54 can be further simplified to the form: 

20 

G(o>) = ( d/2 -^Lm„ *Z n (p) dy 

J-D/2 J0) p - ^» n 



EQ. (57), 

where 2 n (p) is the solution to Sommerfeld's integral. The frequency content of the voltage detected with coil 
assembly 13 will determine the variable G(w). Also, with every value of y and n, the control device 18 may 
30 determine a value of Z n (p) by using well known tables and methods to the various solutions of Sommer- 
feld's integral. Accordingly, the transverse magnetization M T can be determined and an image of the object 
can be generated. 

Another embodiment of the invention is illustrated in Fig. 6. As shown in this figure, the MRI apparatus 
comprises a conventional set of gradient coils 4 and a supply source 8 for generating gradient fields. The 

35 MRI apparatus further comprises the circuitry shown in Fig. 4 of the invention, which generates gradient 
fields by using electric fields. The MRI apparatus of Fig. 6, like the embodiment of Fig. 4, also comprises 
an interface circuit 15, a modulator 17, and an RF amplifier 19. 

The gradient fields produced by the electric fields may be used in the embodiment of Fig. 6 to enhance 
or dephase some of the gradient fields generated by the conventional gradient coils. This may be 

40 advantageous in applications where an image is desired of only a small region of the object. The gradient 
fields produced by the electric fields would combine with the gradient fields produced by the conventional 
gradient coils with the resultant gradient fields affecting only that small region in the object. Thus, a more 
narrow image of the object is produced which reveals only that small desired portion of the object and 
which deletes all other portions of the object, . . 

45 j...TK1b RF coil assembly 13 may also be used to correct for inhomogeneities of the complex dielectric 
constant and also to correct for anisotropy. These errors may result from an inhomogeneous complex 
dielectric constant in the object that shifts the "gradient" horizontally in opposite directions, depending upon 
the polarity of the applied electric field. Thus, the inhomogeneous complex dielectrics distort the magnetic 
field gradients produced by the RF coil assembly 13 thereby rendering the gradient fields non-linear. 

so To illustrate how the RF coil assembly 13 may correct these errors, first consider the expression for the 
phase 4>, which is: 



55 



<t> = yH 0 t + f X f * (yECcosu> D t) dtdr 
J o J 0 

EQ. (58) 
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Defining a first coordinate system T as being equal to r + js, a first detected phase *i is equal to: 

♦* (c> = 5 P E 4o C " dz 



EQ. (59), 



10 

where the variable C (r) is defined accordingly: 

75 



and 

20 



EQ. (60), 



i r D /2 

DJ-D/2 lT) 



25 



EQ. (61). 



After substituting in the expression for C (r) into Equation 59, the expression for the first detected phase *i 
becomes: 



30 

YSin(i) D t 



35 

EQ. (62). 

If we define a second coordinate system W as being equal to u + jv, a second detected phase <b 2 is 
equal to: 

40 

ysin<o D t 



45 



EQ. (63), 



50 



55 



where the variable C (u) is defined accordingly: 



C <u) A C + € (ti) 



and 



EQ. (64), 
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c * -±f C< u) du 

DJ-D/2 {U) 



DJ-D/2 

EQ. (65) . 

Next, after substituting in the expression for C (u) intoEquation 63, the expression for the second 
detected phase <f> 2 becomes: 

70 

ysino) t 



p 

15 



* 2 <c, . 2fS^^<eu./%„,*.> 



EQ, (66). 



If the two coordinate systems T and W are respectively defined by electric fields of opposite polarities, 
which occur when the applied voltages have opposite polarities, then the following variables in the two 
20 coordinate systems are interrelated according to the following expressions: 

T = -W; r = -u; s = -v; E s = -E v EQ. (67). 

With Equation 66, the expression for the second detected phase </>2 may be altered to the form: 

25 

♦,«« . TJ^l£,- V [c(-r).//.,^d<-*)] 

30 

EQ. (68), 

which reduces to the form: 

35 

. Iii^£( £ ,)[cr-/;' e ,. I ,d(-r,] 

40 • 

EQ. (69). 

Since the following equation is true:/ .. 

EQ. (70) , 

the sum of the first detected phase <Jn with the second detected phase 4>2 will be equal to: 



55 
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EQ. (71), 

which is a corrected phase signal that is not influenced by the error signal. Also, the difference between the 
first detected phase 4n and the second detected phase <f> 2 will equal: 

10 

EQ. (72), 

which separates the error signal from the phase encoded signal. 

By subtracting the detected signals when the applied voltages are of opposite polarity to each other, the 

20 errors caused by inhomogeneity and anisotropy in the object cancel. Since the corrected signal therefore 
has less distortion, the image is improved. 

An advantage of the present invention is that the costly and complex gradient coils in the conventional 
MRI apparatus are no longer a necessity since the gradient fields may be produced by the application of a 
time varying voltage V p to the RF coil assembly 13. As a result of the time varying voltage V p , electric fields 

25 are generated in the object and produce time varying magnetic fields that have a uniform curl. Con- 
sequently, the apparatus of the present invention produces linear magnetic field gradients. 

Another advantage of the present invention is that the processing of the detected signals is not limited 
to the two or three dimensional truncated Fourier series or Fourier transformation. Rather, the apparatus of 
the present invention may process the detected signals by projection techniques and may use Bessel 

30 functions or the solution to Sommerfeld's integral in order to solve for the transverse magnetization M T . 

A further advantage of the present invention is that some image artifacts are eliminated. The invention 
does not produce image artifacts caused by the periodic displacement of spins resulting from the 
interaction of their charge and mass with the applied electric field, which is a component of the complex 
permittivity, but is orthogonal to both the z axis and also the gradient field G, since no phase dispersions 

35 are produced unless another gradient is applied along the direction of the applied electric field. Combining 
these gradients as in the embodiment of Fig. 6 allows selective dephasing of some of the spins in the 
imaging volume. 

The image may also be improved with the present invention by detecting phase encoded signals when 
the applied voltages have opposite polarities. By summing together phase encoded signals detected when 
40 the applied voltages have opposite polarities, the error components of the signals cancel each other thereby 
producing a corrected phase encoded signal. The error signal may also be separated from the phase 
encoded signal by subtracting the two phase encoded signals. The corrected phase encoded signal, since it 
no longer contains the error signal, has less distortion which thereby produces an improved image of the 
object. 

45 The foregoing description of the preferred embodiments of the invention has been presented for 
purposes of illustration and description. It is not intended to be exhaustive or to limit the invention to the 
precise forms disclosed. Many modifications and variations are possible in light of the above teaching. 

For example, although the embodiments have been described with reference to a system that uses RF 
pulses, the invention is equally applicable to a continuous wave system. Further, rather than applying a 

50 voltage V p whichhas a constant frequency of o> p , the frequency of the applied voltage V p may vary in order 
to find the resonances of the dispersion of the complex dielectrics of the object and to measure complex 
dielectrics in the object. 

Also, although the electric fields have been used to produce magnetic field gradients for the generation 
of an image, the electric fields may also be used to generate magnetic field gradients which function to 
55 destroy parts of an image. Thus, a more narrow image may be generated which does not reveal portions of 
the object which are not of interest. The technique of making an image more narrow may be used in 
spectroscopy to thereby make it possible to tell the number of spins of a particular species. It may also be 
possible to identify particular species of charged particle mobility with this technique. 

16 

ISDOCID: <EP 06311 51 A1_l_> 



EP 0 631 151 A1 



Moreover, the generation of magnetic field gradients with electric fields has applications with magnetic 
resonance other than in just the imaging of an object. For instance, a magnetic resonance apparatus may 
be used to determine whether fruit is fresh, whether a chocolate bar has been frozen and then thawed, as 
well as a multitude of other possibilities. The invention is equally applicable to these other magnetic 
resonance applications. 

The embodiments were chosen and described in order to best explain the principles of the invention 
and their practical application to thereby enable others skilled in the art to best utilize the invention and 
various embodiments and with various modifications as are suited to the particular use contemplated. It is 
intended that the scope of the invention only be limited by the claims appended hereto. 

Claims 



1, An apparatus for generating magnetic field gradients for use in magnetic resonance imaging of an 
object, comprising: 

is a plurality of conductive wires arranged in parallel to each other and formed in a shape of a 

cylinder; 

a first ring conductor for joining said conductive wires at one end of said conductive wires; 
a second ring conductor for joining said conductive wires at an end opposite said one end of said 
conductive wires; 

20 means for generating a time varying electric field having a first frequency which is lower than a 

Larmor frequency of said object; and 

means for successively applying said time varying electric field between pairs of diametrically 
opposed wires at a second spatial frequency, said second spatial frequency being much less than said 
first frequency; 

25 wherein said time varying electric field produces a linear magnetic field gradient in said object. 

2. An apparatus as set forth in claim 1, wherein said wire conductors are arranged in a bird cage 
configuration that is tuned to said Larmor frequency. 

30 3. An apparatus as set forth in claim 1 or 2, further comprising: 

means for generating RF fields having a frequency centered near said Larmor frequency; and 
" means for applying said RF fields to said conductive wires; 
wherein said RF fields become circularly polarized and rotate about said conductive wires. 

35 4. An apparatus as set forth in claim 1 , 2, or 3, further comprising: 

a first set of filters connected on said first ring conductor between said conductive wires for tuning 
said conductive wires near said Larmor frequency; and 

a second set of filters connected between said applying means and said first ring conductor for 
passing through said first frequency and for filtering out said Larmor frequency. 

40 

5. An apparatus as set forth in any one of claims 1 to 4, wherein said conductive wires detect signals 
generated by said object and said apparatus further comprises means for processing the detected 
signals in order to produce an image of said object. 

45 6. An MRI apparatus for producing an image of an object, comprising: 

means for producing a homogeneous magnetic field along a first axis; 

means for generating a first modulated voltage having a first frequency that is equal to a Larmor 
frequency, said Larmor frequency being specified by said homogeneous magnetic field and by said 
object; 

50 means for generating a second modulated voltage having a second frequency that is less than said 

first frequency; 

means for receiving said first voltage and said second voltage and for generating a first electric 
field and a second electric field respectively having said first frequency and said second frequency; 
control means for selectively applying said first voltage and said second voltage to said receiving 
55 means; and 

means for detecting signals emitted by said object after the application of said first voltage and 
said second voltage; 

wherein said second electric field produces a linear magnetic field gradient in said object and said 
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15 



20 



control means further comprises means for processing the detected signals to generate an image. 

7 ' IS? iKSJJir f ° rth in C ' aim * ^ Said ^ * « Said *- - -ond 

cy.inder ,Ura,ity * C ° ndUCtiVe '" n Para " el t0 each other •* in a shape of a 



30 



35 



- 14 ' iSX'X'SSS ma9netiC fieW 9radiemS ^ ^ h ""»"-» — waging of an object, 

toTninn^H ° f , conductive wires in Parallel to each other in a shape of a cylinder 

^jo.nmg sa.d plurahty of conductive wires with a first ring conductor at one end o, said conductive 

< 5 end I 0 ,' saTd t^X^^ "~ ' *» 31 3 " end said -~ 

freguTy^aidScr^ 9 ^ fie ' d 3 ,irSt WhiCh is ^ «™ a Larmor 

at a s U e ron S H S i V ! ,y ^ ^ Varyi " 9 e ' eCtriC field between P airs <* diametrically opposed wires 

at a second frequency, said second frequency being much less than said first frequenJ 

where.n sa.d t,me varying electric field produces a linear magnetic field gradient ^ said object. 

15 - ,orth in c,aim 14 - ,urther — 9 - °< 

16. A method as set forth in claim 14 or 15, further comprising the steps of- 

generating RF fields at a frequency centered near said Larmor frequency; and 
applying said RF fields to said conductive wires- 

wherein said RF fields become circularly polarized and rotate about said conductive wires. 
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17. A method as set forth in claim 14, 15, or 16, further comprising the steps of: 

connecting a first set of filters on said first ring conductor between said conductive wires for tuning 
said conductive wires near said Larmor frequency; and 

connecting a second set of filters between said applying means and said first ring conductor for 
5 passing through said first frequency and for filtering out said Larmor frequency. 



18. A method as set forth in any one of claims 14 to 17, further comprising the steps of: 

detecting signals generated by said object with said conductive wires; and 
processing the detected signals in order to produce an image of said object. 

10 

19. A method for generating images of an object with magnetic resonance, comprising the steps of: 

producing a homogeneous magnetic field along a first axis; 

generating a first modulated voltage having a first frequency that is equal to a Larmor frequency, 
said Larmor frequency being specified by said homogeneous magnetic field and by said object; 
is generating a second modulated voltage having a second frequency that is less than said first 

frequency; 

receiving said first voltage and said second voltage and generating a first electric field and a 
second electric field respectively having said first frequency and said second frequency; 

selectively applying said first voltage and said second voltage to said receiving means; and 
20 detecting signals emitted by said object after the step of selectively applying said first voltage and 

said second voltage; and 

processing the detected signals and generating an image of said object; 
wherein said second electric field produces a linear magnetic field gradient in said object. 

A method for generating images as set forth in claim 19, further comprising the steps of: 

arranging a plurality of conductive wires in parallel to each other in a shape of a cylinder; 
joining said plurality of conductive wires with a first ring conductor at one end of said conductive 
wires; 

joining said plurality of conductive wires with a second ring conductor at an end opposite said one 
end of said conductive wires; 

wherein said first electric field rotates about said conductive wires at said Larmor frequency and 
said second electric field rotates about said conductive wires at a third frequency that is much less than 
said second frequency. 

35 21. A method as set forth in claim 19 or 20, wherein said step of generating said first and second electric 
fields is performed by an antenna structure that also performs said step of detecting said signals. 

22. A method as set forth in claim 19, 20, or 21, wherein said step of processing uses a projection 
reconstruction technique to produce said image of said object. 

40 

23. A method as set forth in any one of claims 19 to 22, wherein said step of processing employs at least 
one look up table containing Bessel functions when processing said detected signals to produce said 
image of said object. 

45 24. A method as set forth in any one of claims 19 to 23, wherein said step of processing employs at least 
one look up table containing solutions to Sommerfeld's integral when processing said detected signals 
to produce said image of said object. 

A method as set forth in any one of claims 19 to 24, wherein said step of detecting said signals 
comprises the steps of detecting said signals when said second voltage has a first polarity and when 
said second voltage has a polarity opposite said first polarity and said step of processing said signals 
and generating said image comprises the step of adding said signals detected at opposite polarities 
together in order to cancel out an error component in said detected signals. 

55 26. A method as set forth in any one of claims 19 to 25, wherein said step of detecting said signals 
comprises the steps of detecting said signals when said second voltage has a first polarity and when 
said second voltage has a polarity opposite said first polarity and said step of processing said signals 
and generating said image comprises the step of subtracting said signals detected at opposite 
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